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Abstract 

Recent data point in the direction of a cosmological constant dominated universe. We 
investigate the role of supersymmetric QCD with Nf < N c as a possible candidate for 
dynamical cosmological constant ("quintessence"). When Nf > 1, the multiscalar dynamics 
is fully taken into account, showing that a certain degree of flavor symmetry in the initial 
conditions is necessary in order to attain the required late-time "tracking" attractor behavior. 
A possible solution to the problem of long-range interactions of gravitational strength is 
proposed. Finally we show how, as a consequence of the interaction with the Standard Model 
fields, the early cosmological evolution of the scalar fields may be dramatically altered, and 
the approach to the tracking behavior may become much more efficient. 
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1 Introduction 

Indications for an accelerating universe coming from redshift-distance measurements of High- 
Z Supernovae la (SNe la) [|], 0], combined with CMB data || and cluster mass distribution 
0, have recently drawn a great deal of attention on cosmological models with Q m ~ 1/3 
and £l\ ~ 2/3, Q m and Q\ being the fraction densities in matter and cosmological constant, 
respectively. More generally, the role of the cosmological constant in accelerating the universe 
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expansion could be played by any smooth component with negative equation of state Pq/pq = 
wq —0.6 ]|, H, as in the so-called "quintessence" models (QCDM) |J, otherwise known 
as xCDM models @ . 

A natural candidate for quintessence is given by a rolling scalar field Q with potential 
V(Q) and equation of state 

_ Q 2 /2-V(Q) 
WQ Q 2 /2 + V(Q)' 

which - depending on the amount of kinetic energy - could in principle take any value from 
— 1 to +1. The study of scalar field cosmologies has shown || [9| that for certain potentials 
there exist attractor solutions that can be of the "scaling" [jlO], [y|, [T^] or "tracker" |L3|, [L4j 
type; that means that for a wide range of initial conditions the scalar field will rapidly join 
a well defined late time behavior. 

If pq <C pb, where ps is the energy density of the dominant background (radiation or 
matter), the attractor can be studied analytically. 

In the case of an exponential potential, V(Q) ~ exp (— Q) the solution Q ~ hit is, under 
very general conditions, a "scaling" attractor in phase space characterized by Pq/pb ~ 
const [10, 11, 12]. This could potentially solve the so called "cosmic coincidence" problem, 
providing a dynamical explanation for the order of magnitude equality between matter and 
scalar field energy today. Unfortunately, the equation of state for this attractor is wq = wb, 
which cannot explain the acceleration of the universe neither during RD (w ra d = 1/3) nor 
during MD {w m = 0). Moreover, Big Bang nucleosynthesis constrain the field energy density 
to values much smaller than the required ~ 2/3 || |ll], |l2| . 

If instead an inverse power-law potential is considered, V(Q) = M +a Q~ a , with a > 0, 
the attractor solution is Q ~ j^-n/m^ w h ere n _ 3(wq + 1), m = 3(wb + 1); and the equation 
of state turns out to be wq = {wsot — 2)/(a + 2), which is always negative during MD. 
The ratio of the energies is no longer constant but scales as pq/pb ~ a m ~ n thus growing 
during the cosmological evolution, since n < m. pq could then have been safely small during 
nucleosynthesis and have grown lately up to the phenomenologically interesting values. These 
solutions are then good candidates for quintessence and have been denominated "tracker" in 
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The inverse power-law potential does not improve the cosmic coincidence problem with 
respect to the cosmological constant case. Indeed, the scale M has to be fixed from the 
requirement that the scalar energy density today is exactly what is needed. This corresponds 
to choosing the desired tracker path. An important difference exists in this case though. The 
initial conditions for the physical variable pQ can vary between the present critical energy 
density p° cr and the background energy density ps at the time of beginning |l4]] (this range 
can span many tens of orders of magnitude, depending on the initial time), and will anyway 
end on the tracker path before the present epoch, due to the presence of an attractor in phase 
space |l^, On the contrary, in the cosmological constant case, the physical variable p\ 
is fixed once for all at the beginning. This allows us to say that in the quintessence case the 
fine-tuning issue, even if still far from solved, is at least weakened. 

A great effort has recently been devoted to find ways to constrain such models with present 
and future cosmological data in order to distinguish quintessence from A models [15, 16 1. An 
even more ambitious goal is the partial reconstruction of the scalar field potential from 
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measuring the variation of the equation of state with increasing redshift [|17]] . 

On the other hand, the investigation of quintessence models from the particle physics 
point of view is just in a preliminar stage and a realistic model is still missing (see for example 
refs. Jl8| , 19, ^0], ^])- There are two classes of problems; the construction of a field theory 
model with the required scalar potential and the interaction of the quintessence field with the 
standard model (SM) fields [22|. The former problem was already considered by Binetruy 
fj~8|] , who pointed out that scalar inverse power law potentials appear in super symmetric 
QCD theories (SQCD) [^3] with N c colors and Nf < N c flavors. The latter seems the 
toughest. Indeed the quintessence field today has typically a mass of order H ~ 10" 33 eV. 
Then, in general, it would mediate long range interactions of gravitational strength, which 
are phenomenologically unacceptable. 

In this paper we will consider in more detail these problems in the framework of SQCD. 
In our analysis we develope two aspects. First of all, these models have Nf independent 
scalar directions in the vacuum manifold. In ref. [18| all of them were given the same 
initial conditions, so that the dynamics reduced effectively to that of a single scalar field 
with an inverse power law potential of the type considered above. On the other hand, in a 
cosmological setting there is no a priori justification for this assumption, and the fields will 
in general start from different initial conditions. The dynamics is then truly iVj-dimensional 
and it is relevant to know what is the late time behavior of the system, whether there are 
other attractors besides the single-scalar one discussed in ref. [^] , and under what conditions 
is the latter reached by the system. Indeed, we will see that starting with the same initial 
total energy but different initial conditions for the Nf fields may prevent them to reach the 
attractor, so that SQCD cannot be considered as a simple one-scalar model for quintessence 
in these regions of initial conditions phase space. 

The second aspect that we consider is the interaction of the SUSY SU(N C ) model with 
the SM (or supersymmetric SM (SSM)) fields. The quintessence fields are usually assumed 
to be singlets under the SM gauge group and to interact with the rest of the world only 
gravitationally, i.e. via non renormalizable operators suppressed by inverse powers of the 
Planck mass. This is however not enough. In order to prevent long-range interactions of 
gravitational strength it is necessary to assume that the - a priori unknown - couplings 
between the quintessence fields and the SM sector are strongly suppressed today. We do not 
solve this problem, but point out that if a least coupling principle of the type proposed by 
Damour and Polyakov [24] for the superstring dilaton were operative, quintessence models 
could be reconciled with the experimental constraints on the weak equivalence principle and 
on time variation of the SM coupling constants. At the same time, during RD it would be 
quite likely to have SUSY breaking and mass generation for the quintessence fields, with 
masses proportional to H, by the same mechanism discussed by Dine, Randall, and Thomas 
in [25]. If present, these time-dependent SUSY breaking masses would prevent the fields to 
take large values, thus driving the system towards the tracker solution. 

The paper is organized as follows. In Section 2 we review SQCD theories with Nf < 
N c and derive the equation of motion of the scalar degrees of freedom in the cosmological 
background. In Section 3 the tracker solution is discussed. In Section 4 we address the 
problem of the interaction of the superlight fields with the other SM fields and propose a 
possible line of solution. In Section 5 numerical results are presented, with and without the 
interaction with other fields. Finally, in Section 6, we present our conclusions. 
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2 SUSY QCD 



As already noted by Binetruy |18| , supersymmetric QCD theories with N c colors and Nf < N c 
flavors [23] may give an explicit realization of a model for quintessence with an inverse power 
law scalar potential. The remarkable feature of these theories is that the superpotential is 
exactly known non-perturbatively. Moreover, in the range of field values that will be relevant 
for our purposes (see below) quantum corrections to the Kahler potential are under control. 
As a consequence, we can study the scalar potential and the field equations of motion of the 
full quantum theory, without limiting ourselves to the classical approximation. 

The matter content of the theory is given by the chiral superfields Qi and Qi (i = 
1 . . . Nj) transforming according to the N c and N c representations of SU(N C ), respectively. 
In the following, the same symbols will be used for the superfields Qi, Q { , and their scalar 
components. 

Supersymmetry and anomaly-free global symmetries constrain the superpotential to the 
unique exact form 

/ \(3N C -N f ) \ N c —Ns 

where the gauge-invariant matrix superfield Tjj = Qi ■ Qj appears. A is the only mass scale of 
the theory. It is the supersymmetric analogue of Aqcd, the renormalization group invariant 
scale at which the gauge coupling of SU (N c ) becomes non-perturbative. As long as scalar 
field values Qi,Qi 3> A are considered, the theory is in the weak coupling regime and the 
canonical form for the Kahler potential may be assumed. The scalar and fermion matter 
fields have then canonical kinetic terms, and the scalar potential is given by 

_ Nf 1 
ViQuQi) = £ (Vq/ + \Fq} 2 ) + -D a D a (2) 

i=l 

where F Qi = dW/dQi, Fq = dW/dQi, and 

D a = Qlt a Q i -Q i t a 0l. (3) 

The relevant dynamics of the field expectation values takes place along directions in field 
space in which the above D-term vanish, i.e. the perturbatively flat directions (Qi a ) = (£?L)> 
where a = 1 • • • N c is the gauge index. At the non-perturbative level these directions get a 
non vanishing potential from the F-terms in (||), which are zero at any order in perturbation 
theory. Gauge and flavor rotations can be used to diagonalize the (Qi Q ) and put them in the 
form 

in \ _ /nt v _ qA a l < a < N f 

{Wia > ~ Wia> ~ Nf < a < N c 

Along these directions, the scalar potential is given by 

v(qi) = (V(Qi,Qi))-- 




SN C - Nf 1 
N c - N f ' N c - N f ' 
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In the following, we will be interested in the cosmological evolution of the Nf expectation 
values qi, given by 

dV 

(Q l + 3HQ i + —r) = , i = l,...,N f . 
dQl 

In ref. [ |lg| ] the same initial conditions for all the Nf VEV's and their time derivatives were 
chosen. With this very peculiar choice the evolution of the system may be described by a 
single VEV q (which we take real) with equation of motion 



q + 3Hq - g 



A 



2d 



,2g+l 







N c + N 



f 



N c — Nf 



(4) 



thus reproducing exactly the case of a single scalar field in the potential V = A 4+29 ^~ 2g /2 
considered in refs. g, ||, 14]. In this paper we will consider the more general case in which 
different initial conditions are assigned to different VEV's, and the system is described by 
Nf coupled differential equations. Taking for illustration the case Nf = 2, we will have to 
solve the equations 



q\ + 3Hqi - d ■ q% 



q 2 + 3Hq 2 - d ■ q 2 



A 2a 



71 92 J 



A 2a 



\2dN c 



2 + iV c 4 

it 







n 2 

2 + N c % 



- ' "'• ...2 " ' ( 5 ) 

with H 2 = 8"7r/3Mj, (p m + p r + pq), where Mp is the Planck mass, p m i r \ is the matter 
(radiation) energy density, and pq = 2(q 2 + q 2 ) + v(qi, q 2 ) is the total field energy. 



3 The tracker solution 

In analogy with the one-scalar case, we look for power-law solutions of the form 

qtr,i = Ci-t p > , i = 1, • • • , Nf . (6) 

It is straightforward to verify that - when pq <C pp - the only solution of this type is given 
by 



1 - r 



with 



Ci = C 



X 



x l-r A 2(3-r) 

4 m (1 + r) 



1/4 



i = N f 



[1 - r) 2 [12 - m(l + r)] ' 

where we have defined r = Nf/N c {= 1/N C , . . . , 1 — 1/N C ). This solution is characterized by 
an equation of state 

i i 

(7) 



1 + r 



1 -r 



-w B 



2 2 

Eq. (0) can be derived as usual from energy conservation i.e. d (cl 3 pq) = —3 a 2 pQ. 
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Following the same methods employed in ref. || one can show that the above solution 
is the unique stable attractor in the space of solutions of eqs. (H)). Then, even if the q^s 
start with different initial conditions, there is a region in field configuration space such 
that the system evolves towards the equal fields solutions (W), and the late-time behavior is 



indistinguishable from the case considered in ref. [18]. 

The field energy density grows with respect to the matter energy density as 



PQ 3(i+l) 

— ~ a 2 



(8) 



where a is the scale factor of the universe. The scalar field energy will then eventually domi- 
nate and the approximations leading to the scaling solution (§) will drop, so that a numerical 
treatment of the field equations is mandatory in order to describe the phenomenologically 
relevant late-time behavior. 

The scale A can be fixed requiring that the scalar fields are starting to dominate the energy 
density of the universe today and that both have already reached the tracking behavior. The 
two conditions are realized if 



v(qo) 



Pint 



v"(qo) * Hi 



(9) 



where p\ 



crit 



respectively. Eqs. @ imply 



ZMpHfi /8tt and qo are the present critical density and scalar fields VEV 



A 

Mp 



Ml 



3 (l + r)(3 + r) 1 



47r (1 - r) 2 rN c 

3 (l + r)(3 + r) 1 
4^ (1-r) 2 WVc 



l + r 
2(3-r) 



1 P\ 



crit 



2rN c M p 



(10) 
(11) 



Depending on the values for Nf and N c , A and q$/A assume widely different val- 
ues. A takes its lowest possible values in the N c — *■ oo (Nf fixed) limit, where it equals 
4 • 10~ 2 (/i 2 /iV|) 1 / 6 GeV (we have used p° cnt /M P = ( 2 - 5 • 10" 31 /i 1/2 ) 4 ). For fixed N c , instead, 
A increases as Nf goes from 1 to its maximum allowed value, Nf = 1 — N c . For N c 20 and 
Nf close to N c , the scale A exceeds Mp. 

The accuracy of the determination of A given in ( |i~0| ) depends on the present error on the 
measurements of Hq, i.e., typically, SA/A = |5f &Hq/Hq ^0.1. 

In deriving the scalar potential (H) and the field equations (ph we have assumed that the 
system is in the weakly coupled regime, so that the canonical form for the Kahler potential 
may be considered as a good approximation. This condition is satisfied as long as the fields' 
VEVs are much larger than the non-perturbative scale A. From eqs. ( JTo| ) and (O), one can 
compute the ratio between the VEVs today and A, and see that it is greater than unity for 
any Nf as long as N c ^ 20. 

On the other hand, if we want to follow the cosmological evolution of the fields starting 
from an earlier cosmological epoch, we must impose the stronger condition that the q$ have 
been much larger than A throughout the time interval of interest. Taking the tracker solution 
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(lul) clS £1 reference, at the initial redshift Zi n we have, 

z^ 1 -^/ 4 if z - > z 

„m „ „%n „ ^in ^ e <J 11 ^cq 

Qtr _ Qtr,0 Qtr _ Qtr,0 , 

A A q tr ,0 A 3 (r-l)/4 ., 

11 _i ^eq 

where z eq is the redshift at matter-radiation equivalence. For a given r, the condition q™ S> 
A gives an upper bound on Zi n . Taking for instance N c — > oo (Aj fixed) we get Zj n <C 
10 21 (Aj/i) -1 / 3 . In the numerical computations, we will consider initial conditions such that 
the weak coupling regime is always realized. 



4 Interaction with the visible sector 



The superfields Qi and Q { have been taken as singlets under the SM gauge group. Therefore, 
they may interact with the visible sector only gravitationally, i.e. via non-renormalizable 
operators suppressed by inverse powers of the Planck mass, of the form 



d A 9 K j ( 



7i , r 9n ~ "n 



M 2 p n 



M 2 p n 



(12) 



where 4>j represents a generic standard model superfield. From (|10| ) we know that today the 
VEV's qi are typically O(Mp), so there is no reason to limit ourselves to the contributions 
of lowest order in \Q\ 2 /Mp. Rather, we have to consider the full (unknown) functions 



a- 



Q\Qi 

M| 



oo 

£ 

n=0 



ft 



..(q\Qs 



M 2 p n 



and the analogous /?'s for the Qj's. Moreover, the requirement that the scalar fields are on 



the tracking solution today, eqs. (E|) implies that their mass is of order ~ Hq 



10 



-33 



eV. 



The exchange of very light fields gives rise to long-range forces which are constrained by 
tests on the equivalence principle, whereas the time dependence of the VEV's induces a time 



variation of the SM coupling constants [22, 26]. These kind of considerations sets stringent 
bounds on the first derivatives of the /J^'s and (3 l 's today, 



a 



J* 



dlogft* [xf 



dx; 



or — 



dxi 



where Xi = qi/Mp. To give an example, the best bound on the time variation of the fine 
structure constant comes from the Oklo natural reactor. It implies that \6t/a\ < 10~ 15 yr _1 



1 27], leading to the following constraint on the coupling with the kinetic terms of the elec- 
tromagnetic vector superfield V, 



Vi ^-Vi 



a- -, a' ' £ 1(T D i^Mp. 



(ft) 



(13) 
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where {qi) is the average rate of change of qi in the past 2 x 10 9 yr. 

Similar -although generally less stringent™ bounds can be analogously obtained for the 
coupling with the other standard model superfields [^] . Therefore, in order to be phenomeno- 
logically viable, any quintessence model should postulate that all the unknown couplings f3 Jl, s 
and p Jl, s have a common minimum close to the actual value of the (/i'sQ. 

The simplest way to realize this condition would be via the least coupling principle in- 



troduced by Damour and Polyakov for the massless superstring dilaton in ref. 24], where 



a universal coupling between the dilaton and the SM fields was postulated. In the present 
context, we will invoke a similar principle, by postulating that f3^ % = j3 and /J 7 * = /3 for any 
SM field (fij and any flavor i. For simplicity, we will further assume (3 = /3 . 

The decoupling from the visible sector implied by bounds like ( |l3| ) does not necessarily 
mean that the interactions between the quintessence sector and the visible one have always 
been phenomenologically irrelevant. Indeed, during radiation domination the VEVs qi were 
typically <C Mp and then very far from the postulated minimum of the /3's. For such values 
of the qi's the /3's can be approximated as 







M 2 



where the constants (3q and (3\ are not directly constrained by (|T3|). The coupling between 
the (|l4|) and the SM kinetic terms, as in (12), induces a SUSY breaking mass term for the 



scalars of the form 

AL-H 2 ^ (l^| 2 + |g 4 | 2 ) , (15) 

i 

where we have used the fact that during radiation domination (j^j J d^6K^(<f>j, (f)j)^j ~ Prod- 
If present, this term would have a very interesting impact on the cosmological evolution 
of the fields. First of all one should notice that, unlike the usual mass terms with time- 
independent masses considered in [20], a mass m 2 ~ H 2 does not modify the time-dependence 



of the tracking solution, which is still the power-law given in eq. (]6|). Thus, the fine-tuning 
problems induced by the requirement that a soft-supersymmetry breaking mass does not 
spoil the tracking solutions |^0[ are not present here. 

Secondly, since the Q and Q fields do not form an isolated system any more, the equation 
of state of the quintessence fields is not linked to the power-law dependence of the tracking 
solution. Taking into account the interaction with the SM fields, represented by H 2 , we find 
the new equation of state during radiation domination (wb = 1/3) 

W 'V = W Q- ^ 9(1 -r") + 6ft 

where wq was given in eq. @. 

From a phenomenological point of view, the most relevant effect of the presence of mass 
terms like (|l5| ) during radiation domination resides in the fact that they rise the scalar 
potential at large fields values, inducing a (time-dependent) minimum. In absence of such 



i 



An alternative way to suppress long-range interactions, based on an approximate global symmetry, was 



proposed in ref. |22 



s 



terms, if the fields are initially very far from the origin, they are not able to catch up with the 
tracking behavior before the present epoch, and pq always remains much smaller than ps- 
These are the well-known 'undershoot' solutions considered in ref. [14|. Instead, when large 



enough masses ( |l5| ) are present, the fields are quickly driven towards the time-dependent 
minimum and the would-be undershoot solutions reach the tracking behavior in time. 

The same happens for the would-be 'overshoot' solutions, J14[|, in which the fields are 
initially very close to the origin, with an energy density much larger than the tracker one, 
and are subsequently pushed to very large values, from where they will not be able to reach 



the tracking solution before the present epoch. Introducing mass terms like (15) prevents 
the fields to go to very large values, and keeps them closer to the traking solution. 

In other words, the already large region in initial condition phase space leading to late- 
time tracking behavior, will be enlarged to the full phase space. In the next section we will 
discuss numerical results with and without the supersymmetry breaking mass (111 



5 Numerical results 



In this section we illustrate the general results of the previous sections for the particular case 
N f = 2,N C = 6. 




Log(1+z) 



Figure 1: The evolution of the energy densities p of different cosmological components is given as a 
funcion of red-shift. All the energy densities are normalized to the present critical energy density p® r . 
Radiation and matter energy densities are represented by the short-dashed lines, whereas the solid line is 
the energy density of the tracker solution discussed in Section 3. The long-dashed line is the evolution of 
the scalar Geld energy density for a solution that reaches the tracker before the present epoch; while the 
dash-dotted line represents the evolution for a solution that overshoots the tracker to such an extent that it 
has not yet had enough time to re-join the attractor. 



In Fig.l the energy densities vs. redshift are given. We have chosen the same initial 
conditions for the two VEVs, in order to effectively reproduce the one-scalar case of eq. (||), 



9 



already studied in refs. [|, 14]. No interaction with other fields of the type discussed in 
the previous section has been considered. 

We see that, depending on the initial energy density of the scalar fields, the tracker 
solution may (long dashed line) or may not (dash-dotted line) be reached before the present 
epoch. The latter case corresponds to the overshoot solutions discussed in ref. fl4[| , in which 
the initial scalar field energy is larger than ps and the fields are rapidly pushed to very 
large values. The undershoot region, in which the energy density is always lower than the 
tracker one, corresponds to p° CT < p'q < p™. Thus, all together, there are around 35 orders 
of magnitude in p'q at redshift z + 1 = 10 10 for which the tracker solution is reached before 
today. Cleary, the more we go backwards in time, the larger is the allowed initial conditions 
range. 

In Fig. 2 we plot the scalar field equation of state, wq, for the corresponding solutions of 
Fig.l. 
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Fi gUie 2: The cosmological evolution of the equation of state wq = pq/pq for the scalar Geld Q is 
plotted as a function of red-shift. The two cases correspond to the energy densities of the two examples 
in Fig.l. Note that, in the long-dashed curve case, the attractor value of —1/3 for the equation of state 
(corresponding to N c = 6) is joined well before the present epoch and only very recently abandoned, when 
the scalar Geld starts to become the dominating component of our universe and as a consequence its equation 
of state is rapidly driven towards — 1 . 



Next, we explore to which extent the two-field system that we are considering behaves 
as a one scalar model with inverse power-law potential. In fig. 3 we plot solutions with the 
same initial energy density but different ratios between the initial values of the two scalar 
fields. Given any initial energy density such that - for q^/q™ = 1 — the tracker is joined 
before today, there exists always a limiting value for the fields' difference above which the 
attractor is not reached in time. 

The effect of the interaction with other fields discussed in Section 4 is shown in Fig. 4. 
Here, we have included the mass term ( |l5|) during radiation domination with Pi = 0.3 and 
we have followed the same procedure as for Fig.l, searching for undershoot and overshoot 
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solutions. As we see, the range of initial energy densities for the solutions reaching the tracker 
is now enormously enhanced since, as we discussed previously, the fields are now prevented 
from taking too large values. The same conclusion holds even if different initial conditions 
for the two fields are allowed, for the same reason. 

10 32 
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10 20 
10" 
10 14 
10" 
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10 2 
10~ 

2 4 6 8 10 

Log(1+z) 

Fi gUI6 3'. The effect of taking different initial conditions for the fields, at the same initial total field 
energy. Starting with q 1 ™ /qV 1 = 10 the tracker behaviour is not reached today. For comparison we plot 
the solution for q l ™ / q™ = 1 ■ 





Log(1+z) 



Figure 4: The effect of the interaction with other fields, to be compared with Fig. 1 . Adding a term like 
cq. (|15[) with (3\ = 0.3 the would-be overshooting solution (dash-dotted line) reaches the tracker in time. 
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6 Conclusions 



In this paper, we have analyzed the role of SQCD as a possible model for quintessence. Our 
analysis completes the previous one by Binetruy in two respects. 

First, we have taken in full consideration the multi-scalar nature of the model, allowing 
for different initial conditions for the Nf independent scalar VEVs and studying the coupled 
system of Nf equations of motion. Starting with the same initial scalar energy density, but 
different fields' values we have shown that the tracking behavior becomes more difficult to 
reach the larger the difference among the initial conditions for the fields. Thus, an approxi- 
mate flavor symmetry of the initial conditions is needed in order that SQCD may act as an 
effective quintessence model. 

Secondly, we have sketched a possible way out to the common problem of all quintessence 
models considered so far, that is the presence of long-range interactions of gravitational 
strength mediated by the ultra-light scalar fields Our solution is inspired by the Damour- 
Polyakov relaxation mechanism for the superstring dilaton Basically, we postulate that 
all the couplings of the SQCD quark superfields with the SM ones are given by a unique 
function, which has a minimum close to today's values of the scalar fields' VEVs. Since 
all the deviations from Einstein gravity are parametrized by the slope of these couplings 
today, this could make the model phenomenologically safe with respect to limits on the weak 
equivalence principle and on the time dependence of the SM coupling constants. At the 
same time, during radiation domination the coupling with SM fields may have induced a 
SUSY-breaking -time dependent- mass to the scalar fields, with the effect of enhancing the 
initial configuration space leading to a late time tracking behavior. 

A deeper investigation of this scenario requires the implementation of SQCD in a wider 
context, such as superstring theories or theories with large compact extra-dimensions, what 
clearly lies beyond the scope of the present work. 
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